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ABSTRACT 


During the 1955 and 1957 Test Operations at the Nevada Test Site (NTS), masses and 
velocities were determined for more than 20,000 objects, such as glass fragments from 
windows, stones, steel fragments, and spheres, which were energized by blast winds resulting 
from nuclear explosions. Following the field tests, a mathematical model was devised to help 
explain quantitatively the experimental results. This model required certain aerodynamic-drag 
information in regard to the displaced objects. It was the purpose of the study outlined in this 
report to determine the necessary drag properties for the objects by means of drop tests. 

In addition to the objects mentioned above, small laboratory animals, mice, rats, guinea 
pigs, and rabbits, were used in the drop tests. The data obtained from these te^ts were ex¬ 
trapolated to estimate the drag properties for man, and the results compared favorably with 
data from other sources. Also a method was developed to estimate the average drag proper¬ 
ties of man from his total surface area, assuming that every possible orientation of a straight, 
rigid man with respect to the wind was equally likely. 
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NOMENCLATURE 


a = acceleration 
Cd = drag coefficient 

D = sphere diameter or other characteristic dimension 
°F = temperature, degrees Fahrenheit 
Fd = drag force 
F g = force of gravity 
g = acceleration of gravity 
h = height of fall at time t or velocity v 
H = total height of fall 

K = acceleration-coefficient numeric = (p/g) (A/m) (v/D) 


m, M = mass 

m 0 = reference mass 

P 0 = atmospheric pressure, ambient 

Rd = Reynolds number = vD/v 

Ry = a function concerned with error in computed velocity 
R a = a function concerned with error in computed a 
Rq, = an approximation to R a 
R'a - an approximation to R a 
s = frontal area of the object 
t = instantaneous time 
t= thickness of plate-like objects 
t 0 = a reference thickness 
T = total time of fall 

Tnd = total time of fall assuming no drag forces 
AT = T - t N d 
6T = timing error 
v = instantaneous velocity 

V = impact velocity 

V N d = “no-drag" impact velocity 
AV = error in computed impact velocity 
X = time numeric = tg(D/ v) 

AXf= increment of X 

Y = height-of-fall numeric = hg(D/i/) 2 
AY = increment of Y 

Z = velocity numeric = Rd = vD/u 
a = acceleration coefficient = sCj)/m 
a 0 = reference acceleration coefficient 
Aa = error in computed a 

0 T - t/t nd 

6' t = an approximation to 6? 

= an approximation to 6t 

0 V = V / V ND 
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6y = an approximation to 0 V 
v = kinematic viscosity of the air 
p - density of air 
p = density of plate-like objects 
p 0 = a reference density of plate-like objects 
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Chapter 1 

INTRODUCTION 


1.1 OBJECTIVES 

During the 1955 and 1957 Test Operations at the Nevada Test Site (NTS), masses and 
velocities of more than 20,000 objects, 1-5 such as glass fragments from windows, stones, steel 
fragments, and spheres, which were energized by blast winds were experimentally determined, 
and the time-displacement history of an anthropometric dummy simulating man 6 was measured. 
The availability of such a mass of data stimulated an analytical study calculated to arrive at a 
mathematical formulation capable of predicting the translation of objects by blast winds. The 
model 7 designed to predict this translation assumed that certain aerodynamic properties of the 
objects displaced were known. The purpose of this report is to describe the drop-test studies, 
both theoretical and experimental, which were undertaken to determine the aerodynamic 
properties of interest. 

1.2 HISTORICAL BACKGROUND 

When objects dropped in air were first observed to fall at varying rates, it was presumed 
that the difference was inherent in the nature of the materials of which the objects were com¬ 
posed. In the early l§th century, Galileo disproved this theory and established the true cause 
of the observed variations. Small spheres of different substances were dropped at the same 
instant from the Leaning Tower of Pisa, and all hit the ground at nearly the same time. Next, 
Galileo dropped objects composed of the same materials in a variety of shapes and observed 
that they then fell at different rates. Galileo concluded that gravity acted upon all objects with 
the same intensity but that the air resisted their fall in varying degrees, depending on the 
object being dropped. It then followed that all objects would fall at the same rate in a vacuum. 

After the invention of the air pump, Newton verified Galileo’s theory with his famous 
“guinea and feather” experiment. A gold coin, called a “guinea,” and a feather were placed in a 
tube about 5 ft long. When the tube was turned on end, they were observed to fall at different 
rates. When most of the air was removed from the tube with an air pump, no difference in the 
rates of fall of the guinea and the feather was perceivable. 

The force with which the air acts upon a falling object came to be known as a drag force 
since it always resisted the motion of the object. The first experiments to determine the 
nature of this drag force were performed in 1853 by Piobert, Morin, and Didion, 8 who observed 
an object falling in air and recorded its velocity until it reached terrtiinal velocity. Since the 
object at that point was not undergoing acceleration, the drag force equaled the force of gravity. 
They were able to obtain terminal velocities up to about 30 ft/sec. In 1892 Cailletet and 
Colardeau 9 dropped plane surfaces of various shapes from the Eiffel Tower and recorded the 
descent times; they obtained velocities up to 90 ft/sec. In 1905, Eiffel 10 conducted a series of 
experiments that involved attaching the object to be tested to an apparatus that slid down a 
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stretched wire. This apparatus actually recorded both the distance dropped and the drag force 
as a function of time, measured by the vibrations of a tuning fork. Velocities up to 130 ft/sec 
were obtained. 

Today most experiments involving drag measurements are conducted in an artificial air 
stream. 11 This method possesses many advantages over the drop test. In addition to improving 
accuracy, the use of an artificial air stream makes it easier to vary the environment of the 
object being tested. The drop technique is, however, still used to determine air density, 12-15 
to study parachute descents, 16-18 and to determine the drag force on bombs at high subsonic 
speeds. 19 It has also been used to study impact damage to small laboratory animals. 20-23 
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Chapter 2 


EQUATIONS OF MOTION 


2.1 INTRODUC TORY REMARKS 

To predict the motion of an object displaced by the movement of the surrounding air, one 
must know the object’s acceleration coefficient (a). The quantity a has been defined as the 
product of the presenting area and the drag coefficient divided by the mass of the object. 1 The 
same quantity is, of course, involved in the equations of motion of an object in free fall, where 
the object is moving instead of the air. It will be shown in this chapter by development of 
equations of motion that a can be determined experimentally by drop tests if the time of fall, 
height of drop, acceleration of gravity, and air density are known. Analysis of experimental 
data involves the assumptions that (1) no vertical air motion existed, (2) the air density re¬ 
mained constant over the entire drop interval, and (3) the a of the falling object did not change 
during the fall. Measures were taken to reduce air motion during the drop tests (see Chap. 3). 
Since the total drop height was less than 50 ft, the assumption of constant air density is reason¬ 
able. Changes in a may occur during fall by two processes: First, the presenting area may 
change for objects other than spheres because of rotation, and, second, the drag coefficient, a 
function of Reynolds number and thus of velocity of fall, may change. A further discussion of 
these topics will be found in Chap. 4. 


2.2 BASIC EQUATIONS 

A free-falling object is subjected to two forces: a drag force (F^) resulting from the drag 
of the air on the object and the force of gravity (F g ). For the case of free fall in still air, 

Fd = (V2)PV 2 C d s (2.1) 

and 



where F d = drag force 
p = air density 
v = velocity of object 
C D = drag coefficient 
s = frontal area of object 
Fg = force of gravity 
m = mass of object 


g = acceleration of gravity 


( 2 . 2 ) 
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These two forces (Fj and F g ) act in opposition. By Newton’s second law of motion, their 
algebraic sum is equal to the product of the mass and the acceleration (a) of the object, or 

ma = mg-(V 2 )pv 2 C D s 

which can be written in the form: 


a = g-(V 2 )pv 2 


(2.3) 


The quantity sC D /m is the acceleration coefficient (a). Its units are area/mass (ft 2 /lb is 
used throughout this report). Thus 


sCd 


(2.4) 


Alpha has a special significance since it combines the significant parameters of the object 
which influence the drag force and consequently the acceleration of the object, i.e., if two 
objects with the same a were dropped at the same time under identical conditions, their accel¬ 
erations and velocities would remain the same throughout the fall. 

Therefore Eq. 2.3 can now be written as 


= g(i - 


.2 pa 


(2.5) 


By definition, v and a can be written in differential form: 


where t = time and h - height. 
By combining Eqs. 2.6 and 2.7, 


( 2 . 6 ) 


(2.7) 


dh = 


v dv 


( 2 . 8 ) 


and substituting a from Eq. 2.5 into Eq. 2.8, 


dh = 


(>- 


v dv 
2 POl' 


(2.9) 


Now, if a is assumed to be constant throughout the low-velocity range v* < v < v 2 , Eq. 2.9 
can be integrated (see Appendix, Sec. A.1.1) over the range h t < h < h 2 and t t < t < t 2 to give the 
following equations: 

h 2 - hi = (pa)” 1 In {[1 - (pa/2g)^ 2 ] [1 + v x (pa/2g) /z ]} 

- (pa)" 1 In {[1 - v 2 (pa/2g)^ 2 ] [1 + v 2 (pa/2g)^ 2 ]} (2.10) 
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and 


t 2 - tj = (2gpa) y2 In {[1 - vn (pa/2g) 1/2 ] [1 + v 2 (pa/2g) 1/2 ]} 

- (2gpa) - ^ 2 In {[1 - v 2 (pa/2g)^] [1 + Vj (pa/2g)‘ /2 ]} (2.11) 

Before these equations were used for numerical computations, they were transformed to a 
more convenient nondimensional form, which is explained in the Appendix (Sec. A. 1.2). Fig¬ 
ure 2.1 illustrates how the drag coefficient of a sphere varies with Reynolds number. 2 The 
Reynolds number (R d ) is a dimensionless quantity defined by 

vD 

R d=— ( 2 . 12 ) 

where R d = Reynolds number 

v = velocity of the sphere 
D - diameter of the sphere 
v - kinematic viscosity of the air 

Figure 2.2 shows the Reynolds number as a function of velocity and characteristic dimen¬ 
sion (diameter for a sphere, length of side for a square plate, etc.). The characteristic dimen¬ 
sion is read directly in inches; whereas, velocity is read directly only if the kinematic viscosity 
of the air is 0.000179 ft 2 /sec. This value was chosen because it represents the average kine¬ 
matic viscosity at 5300 ft above sea level, 3 the approximate altitude of the southeast section of 
Albuquerque, N. Mex., where the experiments were performed. If the kinematic viscosity has 
any other value, the velocity variable must be modified by the appropriate factor, as indicated 
in Fig. 2.2, to convert it to its true velocity. If the velocity, diameter of a sphere, and the kine¬ 
matic viscosity of the surrounding medium are known, Fig. 2.2 can be used to find the Reynolds 
number. Then for that Reynolds number, the corresponding drag coefficient can be determined 
with the aid of Fig. 2.1. 

2.3 EQUATIONS FOR CONSTANT a 

If a is assumed to be constant throughout a fall; if the initial conditions are v A = hj = = 0; 

and if v 2 = V, h 2 = H, and t 2 = T (see Appendix, Sec. A.1.1), Eqs. A.3 and A.5 become 

V = (2gH) 1/2 [(1 - e~P aH )/panf 2 (2.13) 

and 

T = [(2H/g) V!! (pa H) y2 /2] + (2H/g)‘ /2 (paH)- 1/2 In [1 + (1 - e -p “ H )’ /2 ] (2.14) 

An alternative way of writing Eq. 2.14 is 

T = (2H/g) y2 (paH)- V2 cosh -1 e paH/2 (2.15) 

An explicit solution of Eq. 2.13 for a in terms of V and H was not possible; however, for a 
fixed value of H, V can be plotted as a function of a. Such a plot, illustrated in Fig. 2.3, was 
found to be useful in a previous study of the impact loading of small anesthetized laboratory 
animals 4 in which impact velocity was determined by means of high-speed photography. The 
experimental data shown on this chart (Fig. 2.3), however, were determined in the present 
study (Sec. 4.5). 

Acceleration as a function of height of fall can be found by substituting V from Eq. 2.13 
into Eq. 2.5 to give 

a = g e~ p “ H (2.16) 
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2.4 EQUATIONS WITHOUT DRAG * 

The equations for free fall with no-drag force and initial conditions V = H = 0 when T = 0 

are 

a ND = g (2.17) 

V ND = S t ND (2.18) 


and 

H = g (T nd ) 2 /2 

Equation 2.19 can be solved for T to give 
T ND = (2H/g) 1/2 

Substituting T from Eq. 2.20 into Eq. 2.18 gives 
V ND = (2gH) V2 


(2.19) 

( 2 . 20 ) 

( 2 . 21 ) 


2.5 COMPARISON OF EQUATIONS OF MOTION WITH AND WITHOUT DRAG 
A summary of the preceding equations is given here: 

For any initial conditions: 

(a) With drag (for constant a) 

a. = g (i - v 2 

(b) No drag 

a ND = g 

For initial conditions v = h - 0 when t = 0: 

(a) With drag (for constant a?) 

V = (2gH) V2 [(1 - e-P“ H )/paH] y2 
a = g e"P“ H 

T = [:(2H/g)‘ /2 (paH)*/2] + (2H/g) V2 (paH)-' /2 In [1 + (1 - e-P“ H ) % ] 

(b) No drag 

V ND = (2gH) 1/2 

a ND = g 
T nd = (2H/g) V2 


(2.5) 

(2.17) 

(2.13) 
(2.16) 

(2.14) 

( 2 . 21 ) 

(2.17) 

( 2 . 20 ) 


*The subscript ND is used to indicate the no-drag value of the parameter; i.e., Tnd is 
the time an object would take to fall, assuming the object did not experience drag force. 
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For the initial conditions v = h = 0 when t = 0, the quantities with drag are arranged so 
that they differ from those without drag by a factor that is a function of the dimensionless 
quantity pa H. 

Figure 2.4 shows these results in a different form. Theta (0) is defined as the ratio of a 
variable (V, T, or a)-with drag to the same variable without drag and is a dimensionless 
quantity expressed as a function of the paH. For example, 9y is the final or impact velocity 
divided by the velocity for a drop with no drag. The final or impact velocity can be computed 
from Fig. 2.4 and measured values of H, g, and T (assuming a, g, and p are constant). This 
calculation can be made by computing Tnd and 9 x and then using the 9^ curve to determine 
pa H. For that value of paK, 9y is then read and multiplied by the computed value of Vnd to 
determine the impact velocity. If the values for p and a are known and H is measured, 9 t and 
9 y can be read directly from Fig. 2.4. 

The dotted curves in Fig. 2.4, which are approximations to the solid curves, are good over 
certain ranges of paH and are useful to simplify computations. They were computed from the 
following equations: 

d' T = [(paH) 1/2 /2] + (pa H)~ % In 2 (2.22) 


and 


9(j = (paH) - ^ (2.23) 

These approximate relations depend upon the fact that e _ ‘ v: ' " — 0 very fast as paH -*■ 00 . Hence 
zero was substituted for e“in the equations for and 9y (Fig. 2.4) to arrive at the ap¬ 
proximations. Using Eq. A.16 from the Appendix (Sec. A.2) and the definition 9 v = V/Vnd , one 
can write 

V/Vt = 0V (2gH) ,/2 (2g/pa) - ^ 2 

= 9 V (paH) Vz = dy/e’y (2.24) 

where Vt is the terminal or limiting velocity as pa H —* °°. Thus the convergence of 0y to 0y- 
indicates V Vt or the approach of the velocity of an object to terminal velocity as the 
height of fall increases. Another approximation to 

6>" = 1 + (paH/12) (2.25) 

for low values of paH was found by expanding 9> p in a Maclaurin’s series as follows: 

6>T(paH) = 0 T lo + de T /d(paH) | 0 (paH) + d 2 @ T /d(paH) 2 |„ (paH) 2 /2!+ ... 

= 1 + (paH/12) + (paH) 2 /480 + ... (2.26) 

The approximation shown by Eq. 2.25 is then the first two terms of the Maclaurin’s series for 
0 T (paH). 

Although impact velocity can be computed from height of fall with Fig. 2.4, as described 
in the preceding paragraphs, Fig. 2.5 represents a simpler way to accomplish the same thing. 

It is especially useful if g = 32.174 ft/sec 2 , in which case, V, T, and H can be read directly. 

This value for g was chosen because it represents the acceleration of gravity 5 at sea level. If 
the acceleration of gravity has any other value, the velocity and time variables must be modi¬ 
fied by the appropriate factors, as indicated in Fig. 2.5, to convert them to actual velocity and 
time. 

2.6 ESTIMATIONS OF ERRORS 

To estimate the percentage of error in computed impact velocity (V) introduced by a 
timing error (6T), the following equation was used: 
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(2.27) 


100AV/V = 100(AV/AT)(1/V) at = 100(a/V) 6T 

Substituting for a and V from Eqs. 2.16, 2.13, and 2.20 in Eq. 2.27 gives 

100AV/V = [(paH)^e -paH (1 - e- paH f 1/2 ] 1006T/T ND 

= R V • 1006T/T ND (2.28) 

where Ry is the dimensionless quantity in brackets. Figure 2.6 is a plot of R v vs. paH. To 
find the estimated percentage of error in computed impact velocity due to error in measured 
time of fall, determine Ry from Fig. 2.6 and use this value in Eq. 2.28. It is interesting to 
note that 100AV/V <1006T/T ND for all values of pa H; indeed, by substituting 

T nd = T{[(paH)^/2] + (paH)- V2 In [1 + (1 - e-P« H ) V2 ] -1 } 

from Eqs. 2.14 and 2.20 into Eq. 2.28, it can be shown that AV/V < 6T/T for all values of pa H. 
Thus impact velocity can be determined with greater accuracy, percentagewise, than the 
measurement of time of drop from which it is computed. 

A similar procedure was used to estimate the percentage of error in computed a intro¬ 
duced by a timing error 6T. 

lOOAa/a = 100(Aq:/AT)(1/q')AT = 100 (^ «) 5T ( 2 * 29 ) 

Differentiate Eq. 2.14 with respect to a and substitute in Eq. 2.29 to obtain 

100Aa/a! = 2{[(paH) 1/2 /2] -(paH) -I/2 In [1 + (1 - e - ' 3 ® 11 )^ + (paH)' A e-P««[(l - e-P“H) 

+ (1 - e - paH ) V ' 2 ]- 1 }- 1 1006T/T ND 

= R a 1006T/T ND (2.30) 

where the dimensionless quantity in braces is given the symbol Rq. 

Figure 2.7 is a plot of R a vs. pa H. To find the estimated percentage of error in com¬ 
puted a resulting from error in measured time of fall, determine R a from Fig. 2.7 and use 
this value in Eq. 2.30. The two dashed curves are reasonable approximations to the solid 
curve over certain ranges of pa H. The R^ approximation was formed in the same way the 6j 
and Qy approximations were formed (Eqs. 2.22 and 2.23). The R^ approximation for low 
values of pa H was obtained by differentiating Eq. 2.25 and substituting it into Eq. 2.29 to obtain 

R’ a = 12/paH ( 2 - 31 ) 

Equation 2.31 was used to approximate the timing accuracy needed for the experiment, 
where p « 0.06 lb/ft 3 , H * 50 ft, g w 32 ft/sec 2 , and a w 0.03 ft 2 /lb. Since the timing require¬ 
ment is greatest for low values of a , the value used in this report was the lowest one antici¬ 
pated for the objects to be dropped. From Eq. 2.31, Rq, « 12/pQ!H = 12/0.09. If this value of 
Rq, is used for the quantity in braces in Eq. 2.30, 

lOOAa/a = (12/0.09)1006T[2(50)/32]“tt « 75006T 

Thus, if 5T = 0.2 msec, then lOOAa/a « 1.5 per cent. Therefore to limit the error in com¬ 
puted a to about 1 per cent, the timing mechanism would need to determine the drop time to 
the nearest 0.1 msec. 
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Fig. 2.1— Drag coefficient of spheres vs. Reynolds number. Over a large range of Reynolds numbers, Cd is nearly con¬ 
stant at 0.47. (Data from Ref. 2.) 
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Fig. 2.2— Reynolds number vs. velocity, characteristic dimensions, and kinematic viscosity. Characteristic dimension 
is the diameter of a sphere, length of side of a square plate, etc. 



Fig. 2.3 — Impact velocity vs. height of fall and a. Included in the figure are the animal data from the 
present report showing mean values (represented by o) and standard deviations (represented by ** - *) 
of the experimentally determined a. 





















Fig. 2.5—Impact velocity vs. time of fall computed for several drop heights. Velocity and time can be 
evaluated for any acceleration of gravity. H lines are included for orientation. 








. Find Ry from curve and percentage of error in V from equation shown. 




Fig. 2.7 — Method for estimating error in computed a resulting from error in measured time of fall 
Find R r from curve and percentage of error in a from equation shown. 






Chapter 3 


EXPERIMENTAL PROCEDURE 


The experimental system used to measure the time required for various objects to fall a 
measured distance is shown diagrammatically in Fig. 3.1. The object to be dropped was placed 
near the end of a bar that was designed to swing about its mid-point (see drop device in upper 
right of Fig. 3.1). The drops were initiated manually by tripping the catch shown just left of 
the microswitch. This action allowed a helical spring to rotate the bar about its mid-point. As 
the bar began to move, the test object was left free to fall; and, at the same time, the micro¬ 
switch was activated, starting an electronic timer (Hewlett Packard Electronic Counter, model 
522-B). At the termination of the fall, the test object struck an aluminum plate, creating a 
noise, which was detected by a contact microphone placed on the underside of the impact plate. 
The signal from the microphone then stopped the electronic timer. 

Figure 3.2 is a picture of the drop apparatus and a close-up view of the release mecha¬ 
nism. The dropping bar, supported on both sides by vertical members, was made thin ( 3 / 8 in.) 
to minimize air disturbance during release of the test object. Plates of various sizes were 
made to be attached to the right end of the dropping bar, whose total length was 32 in., so that 
objects as large as a rabbit could be conveniently dropped. A helical spring, hidden from view 
by the vertical support, was designed to ensure that the downward acceleration of the right 
side of the bar would be greater than the acceleration of gravity. After the bar moved approxi¬ 
mately 90°, its energy was partially dissipated by impact with a small pad of Ensolite (a 
product of the U. S. Rubber Co.) placed between the upright supports; thus rebound was less¬ 
ened, preventing interference of the bar with the falling object. 

It was necessary to pay particular attention to the design and construction of the release 
mechanism. If any motion had been imparted to the dropping bar during the process of release, 
then the times of drop would not have been reproducible. Therefore, the sliding parts of the 
latch were made of blocks of tool steel, precision ground and clamped at the ends of the drop¬ 
ping bar; and the vertical bar, shown on the left side of Fig. 3.2, was designed to pivot about a 
point near the base of the apparatus. Even after taking these precautions in design, it was 
found that more reproducible times of drop could be obtained if the catch were released slowly 
by means of the knurled nut shown at the left in Fig. 3.2, thus minimizing motion induced in 
the dropping bar. 

Tests were conducted in an elevator shaft in which objects were dropped 48.372 ft. Pre¬ 
cautions were taken to minimize vertical air currents by keeping all access doors to the shaft 
closed. 

The impact plate shown schematically in Fig. 3.1 consisted of a V 4 -in.-thick 2-ft 2 alumi¬ 
num plate. The plate was fastened to a table-like support about 6 in. high made of 3 / 4 -in.-thick 
plywood in which a hole was made in the center to accommodate the contact microphone. It is 
pertinent to note that the impact plate made a relatively small target for irregular objects 
dropping 48 ft. Occasional “misses”*were experienced, especially by the smaller asymmetri¬ 
cal objects. Plate-like objects dropped flat also tended to veer away from the impact plate if 
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a small degree of rotation occurred. Surprisingly enough, the release operator could usually 
observe the object's fall well enough to detect the presence or absence of rotation. Thus the 
data included in this report are for objects that had little or no rotation while falling. 



MICROPHONE 


Yi g. 3.1 Block diagram of the drop-test equipment showing release apparatus and starting 

and stopping circuits. 
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Fig. 3.2 — Release apparatus used in the drop tests. The helical spring designed to accelerate the motion of the dropping 
bar (shown in horizontal position) is hidden from view behind the vertical support. The inset is an enlargement of the 
release mechanism (see text). 










Chapter 4 


ANALYSIS OF DATA 


4.1 CONSTANTS FOR THE DROP TESTS 


The constants for the drop tests were 
g = 32.126 ft/sec 2 
H = 48.372 ft 
p = 0.060437 lb/ft 3 
v = 0.0002071 ft 2 /sec 

The value of g was obtained from the Coast and Geodetic Survey, Department of Com¬ 
merce. 1 H was measured directly; whereas, p and v were computed by measuring the air 
pressure and temperature and substituting them in the following derived equations: 2 




459.69 + t(°F) 



2.7803 (10)- 7 [456.69 + t(°F)] 


(4.1) 

(4.2) 


Since the drop tests were conducted over a period of three days, the above values of p 
and v are averages; however, the deviations about these average values from day to day were 
small, the range in both cases being less than 1 per cent. 


4.2 TIME CORRECTIONS 

Steel spheres, 7 /i6 and Vie in * in diameter, were dropped during each of the three days when 
drop tests were conducted. Although the drop times for a particular sphere on a given day 
were quite consistent (low standard error), the mean drop time varied significantly from day 
to day. Since the air density remained nearly constant, the differences were probably due to 
minor adjustments that had been made on the microswitch and the release apparatus between 
daily tests. 

The circumstance stated above pointed to the need for accurately computed drop times 
for the two “calibration” spheres in order to determine the timing deviations inherent in the 
experimental equipment. This was accomplished by the stepwise integration of Eqs. 2.10 and 
2.11, the height of drop and time of drop being computed for small steps in velocity for which 
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C D (and thus a) could be considered to have a constant value (see Fig. 2.1). Thus the following 
values were computed for time above no-drag time (AT): For a 3 / 16 -in.-diameter steel sphere, 
AT = 35.05 msec; and for a 7 / 16 -in.-diameter steel sphere, AT = 14.91 msec. 

These values, when compared to the measured times of drop, indicated that the timing 
error for each day was essentially the same for the two spheres. The average corrections to 
be added to the measured times of fall were calculated to be: first day, +3.45 msec; second 
day, -0.14 msec; and third day, -1.65 msec. 

In addition to these corrections, another correction was needed for the tests when heavy 
objects were dropped on a pad of Ensolite. This 1-in.-thick pad was used on the third day to 
help protect the contact microphone. Several spheres were dropped both with and without the 
pad. The differences in times of fall produced by the pad were plotted against the uncorrected 
times of fall with the pad. The points indicated a linear relation. This correction was then 
combined with the constant correction previously determined for the third day to give a total 
time correction (Fig. 4.1). 

4.3 DATA FOR SPHERES, STONES, AND MILITARY DEBRIS 

4.3.1 Introductory Remarks 

It was implied in Sec. 4.2 that to determine an accurate time of drop relation for spheres, 
it was necessary to allow for a varying a dependent on the Reynolds number, which changes 
throughout the fall. If this were true for irregular objects, then the drop-test technique would 
be doomed to failure since the relation between drag coefficient and Reynolds number is not 
known for objects of interest to this study. In this section it will be shown that the dependence 
of a upon Reynolds number for irregular stones is evidently quite different from that for 
spheres. This difference makes it reasonable to assume for experimental purposes that a for 
stones can be considered to be constant for the entire drop interval. 

4.3.2 Sphere Data 

First, it will be necessary to examine the test data for the spheres. Figure 4.2 shows a 
plot of “theoretical” a (based on C D - 0 .47) for spheres vs. measured time of drop in excess 
of the no-drag time. The curve shown was computed assuming a constant a by using Eq. 2.14. 
For most of the spheres, the times measured were either too low or too high when compared 
to the constant a curve. An examination of the lower curve in Fig. 4.3 will help explain those 
deviations. Marked on this chart are the Reynolds numbers computed for the final or impact 
velocities for various spheres. The position of the Vi 6 -in.-diameter steel sphere, for instance, 
indicates that high drag coefficients existed for a large percentage of the fall. The effect of 
this was that the measured time of drop was longer than that expected assuming a constant 
drag coefficient of 0.47 (Fig. 4.2). Times for the V 8 -in. -diameter steel spheres, on the other 
hand, corresponded reasonably well with the constant a assumption. Undoubtedly the reason 
for this was that the effect of the initially high drag coefficients were compensated for by co¬ 
efficients lower than 0.47 near the termination of the drop 3 (Fig. 4.3). For spheres still larger 
in diameter (e.g., 3 / 4 -in.-diameter aluminum spheres), the drop times were less than expected 
since a larger percentage of the drop-interval involved low drag coefficients (Fig. 4.3). The 
sphere drop involving the highest Reynolds numbers was that of a 3.34-in.-diameter croquet 
ball whose weight had been increased with a brass plug. In this instance a large portion of the 
drop entailed drag coefficients whose values were near 0.47. Thus the measured drop time 
corresponded closely to that expected for a constant coefficient of this value. It is interesting 
to note that the “effective” drag coefficient for the entire drop interval can be approximated in 
each individual instance by the coefficient corresponding to one-half the Reynolds number. 

4.3.3 Stones Compared to Spheres 

To study the differences in the Reynolds number-drag relation between spheres and 
irregular objects, it was necessary to analyze the two groups of data similarly. This involved 
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computing effective a for each sphere drop, assuming that the drag coefficient was constant for 
the entire drop (Eq. 2.14). The a values so computed were then plotted as a function of sphere 
mass, as illustrated in Fig. 4.4. If the drag coefficient is assumed to be constant, it can be 
shown that a for spheres of a constant density varies inversely as mass to the / 3 power. Thus, 
if the drag coefficient had the same constant value for all sizes of like-density spheres, the 
least-squares regression lines shown in Fig. 4.4 would have slopes of - Deviations from 
this value indicate that the smaller spheres had effective drag coefficients somewhat greater 
than those for the larger spheres. 

Statistical data pertinent to the slopes of the least-squares regression lines for spheres 
and stones are recorded in Table 4.1. In addition to the slopes, the standard deviations of the 
slopes were computed \ and from these data the probability (pr) that the measured slopes could 
have been deviations from the “true” slope of - */ 3 was computed. Although the probabilities for 
the individual cases varied considerably, that for a combination of all spheres proved to be only 
0.034. Stated in physical terms, the Reynolds number-drag relation for spheres is such that, 
for the purposes of the present study, the effective drag coefficient varies from one size sphere 

to another. 

Table 4.1 also presents the results of a similar statistical treatment for six types of ir¬ 
regular stones. With one exception, pumice, the probabilities that the measured slopes could 
have been deviations from the “true” slope of were quite high. This finding indicated that 
the computed slopes are likely to be random variations from the true slope. The variation of 
the pumice data from the — V 3 slope was probably due to the absorption of moisture from the 
air in the outer layer. If this layer were of uniform thickness, the smaller stones would suffer 
a larger relative increase in weight than would the larger stones. 

4.3.4 Stone Data 

The reasoning presented above seemed to justify a constant-a assumption. The a values 
for individual objects were, therefore, determined graphically from a curve of a vs. AT similar 
to that shown in Fig. 4.2. 

Figure 4.5 is a plot of experimentally determined a for 41 pieces of painted gravel as a 
function of mass. These stones were taken from a sample of stones that was marked for identi¬ 
fication and used in secondary-missile experiments in Nevada during the Operation Plumbbob 
weapons tests. 4 Since the stones actually used in these experiments had masses less than 1500 
mg, a separate analysis was made for these. In addition, another analysis was made for all 
stones, assuming a slope of the least--squares regression line of —V 3 (Sec. 4.3.3). The results 
of both analyses are shown in Fig. 4.5. 

Drop tests were made with limestbne fragments so that there could be a comparison of the 
results with stones whose densities were less uniform. The slope of the least-squares regres¬ 
sion line determined for these data (Fig. 4.6) is reasonably near the expected value of -% 

The data shown in Figs. 4.7, 4.8, and 4.9 are for stones obtained at three different loca¬ 
tions (10P4, 1S1A, and 4.3 GTS) at NTS during Operation Plumbbob. 4 An analysis was made of 
the combined data for stones at the 10P4 and 1S1A locations because the individual analyses 
were remarkably similar. 

Figure 4.10 contains the data for the pumice fragments that were discussed in Sec. 4.3.3. 
The least-squares regression lines for all the stone data are displayed in Fig. 4.11. An exami¬ 
nation of this figure points to the similarity between pumice and stones at the 4.3 GTS location, 
both having higher a for given stone masses than the other samples. 

4.3.5 Military-debris Data 

The military debris used in the drop tests consisted of steel shrapnel and a few other ir¬ 
regular pieces of steel, all of which had been used and recovered as secondary missiles in the 
Operation Plumbbob study. 4 The a-mass data for these steel fragments are shown in Fig. 

4.12. For comparison purposes, the theoretical cn-mass relation for steel spheres, assuming 
C D = 0.47, is shown on the same chart. It is apparent and reasonable that, for a given mass, a 
fragment of steel has a much larger ct than does a steel sphere. The slope of the least-squares 
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regression line is less than that for the theoretical line for spheres; however, it is not as flat 
as a line for plates (Sec. 4.4). 


4.4 DATA FOR PLATES, CUBES, AND GLASS FRAGMENTS 

4.4.1 Introductory Remarks 

Drag properties of window- and plate-glass fragments were of interest to this study be¬ 
cause they constitute biological hazards when hurled through the air by the action of a blast 
wave on a structure containing windows. 4 Regular plates and cubes were studied only for the 
purpose of better interpretation of the glass-fragment data obtained in full-scale nuclear tests. 

4.4.2 Regular Plates and Cubes 

It should be pointed out that the plates referred to in this study are really “thick” plates, 
i.e., no attempt was made to sharpen the edges in order to achieve the “thin” plate that is 
usually used in wind-tunnel studies. For this reason the Reynolds number-drag relation 
shown in the upper part of Fig. 4.3 does not strictly apply to the experimental data reported 
here. However, it is useful to note in Fig. 4.3 that within the region of interest the variation in 
drag coefficient is less than 5 per cent (from 1.13 at R^ = 2 x 10 4 to 1.17 at Rd = 8 x 10 4 ). 

Figure 4.13 presents drag coefficients for regular plates in normal flow determined from 
drop-test data, assuming that a was constant throughout the fall. If the plates used had been 
thin plates, those with the smallest area could have been expected to have the highest drag 
coefficients (Fig. 4.3). Since the reverse was true, it can be concluded that the thickness effect 
must outweigh the effect of Reynolds number. It is interesting to note that no difference in C D 
could be detected between plates that were square, circular, or triangular (equilateral). 

Similar data are presented in Fig. 4.14 for cubes made of various materials. Although a 
considerable change in terminal Reynolds number from the smallest to the largest cubes was 
noticed, no significant trend in the drag coefficient data was apparent. 

With information obtained from the same drop tests for cubes, a was plotted as a function 
of the length of the cube on logarithmic paper (Fig. 4.15). If drag coefficient were constant, 
then the following relation can be stated for cubes made of material of constant density: 
a ~ (s/M) ~ (length of side)^ 1 . Thus a comparison of the data points in Fig. 4.15 to the lines 
drawn with a slope of —1 indicate the degree to which drag coefficient is constant for the 
different types of cubes. 

4.4.3 Data for Glass Fragments 

Two types of glass fragments were used in this study: those originating from window 
panes (0.125 in. thick) and those from plate-glass panes (0.225 in. thick). Figure 4.16 illus¬ 
trates the random shapes of typical window-glass fragments. The number appearing with each 
fragment indicates its mass in grams. 

Drop times obtained for 380 window-glass fragments were used to determine individual a 
by the method previously described, assuming constant a throughout the fall. When a was 
plotted as a function of mass on logarithmic paper, it was apparent that the points tended to 
form not one, but several line segments. The following procedure was used to determine the 
mass boundaries for each line segment. The data were divided into small mass groups, and 
the mean a and standard deviation in a were determined for each group. These data were 
plotted as illustrated in Fig. 4.17. The mass limits for each line segment were then estimated, 
and these limits used to determine the least-squares regression equations for each segment. 
These segments are shown in Fig. 4.18, along with the data points used in their determination. 
Because the smaller fragments were approximately spherical, the regression line for them 
was given a theoretical slope of — % (Sec. 4.3.3). Another group of window glass was dropped 
edgewise (shown as triangles in Fig. 4.18). The regression line for these data passed near the 
intersection of the first two segments. For the sake of simplicity, another regression line 
was computed for the data, with the requirement that it pass through the intersection of the 
first two segments. The standard error of estimate was increased very little by this procedure. 
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Alpha data for plate glass (0.225 in. thick) dropped flat are presented in Fig. 4.19. It is 
interesting to note that a larger scatter is shown in the data for small fragments than for large. 
This, no doubt, can be attributed to the greater mechanical stability during fall of large plates 
than of small ones. This scatter also was observed at the time of the experiments. The same 
distribution of scatter can be detected in the data for window glass dropped flat as illustrated 
in Fig. 4.18. 

A summary of all glass-fragment data is shown in Fig. 4.20. The data represented by line 
segment 1 for the smaller fragments can be used for either plate or window glass in any 
orientation since the fragments were small enough that their flat dimensions were not ap¬ 
preciably different from their thickness dimensions. Thus for plate glass line segment 1 was 
extended to higher masses until it intersected line segment 6 for plate-glass fragments dropped 
flat. This extension was possible because the densities of the window and plate glass were the 
same. 

4.4.4 Analysis for Thick Plates of Random Shape 

The purpose of this section is to derive relations that make it possible to determine drag 
coefficient (and thus a) for plates of any thickness and presenting area within the limits of the 
experimental data. The results apply to plates of random shapes such as those which have 
been described (see Fig. 4.16). The assumption is made, and verified by the experimental 
data, that plates with the same characteristic length-to-thickness ratio have the same drag 
coefficient, at least for the experimental conditions being considered here. 

The characteristic length of a plate of random shape was taken to be the length of a side of 
a square plate of an equivalent area, i.e., the square root of the area. 

It has been shown (see Figs. 4.19 and 4.20) for wide ranges of masses and for plates of 
constant thickness and density that the relation between the acceleration coefficient and mass 
can be expressed by an equation of the form 

log a; 0 = J + B log m 0 (4.3) 

where J and B are constants. With this expression, the relation between Cq and s^/t, where s 
represents the area of a plate in normal flow and t represents its thickness, can be determined. 

The subscript zero is used to designate the parameters of the plates used in this study. 

The same terms without subscripts represent quantities that can vary between limits deter¬ 
mined by the experimental data. Thus a , the acceleration coefficient, may be defined as = 
s 0 CD/m 0 and a = sCd/hi (Sec. 2.2), where the drag coefficients (Cd) are identical for the con¬ 
dition that s 0 1/2 /to = s Vt. For each value of s^/t and, therefore, for each value of C D , it can be 
stated that s = k 4 t 2 p, s 0 = k^Jpo, m = k 2 t 3 p, and m 0 = k 2 t 3 p 0 , where kj and k 2 are constants, t is 
the thickness of the plate, and p is the density of the material from which the plate is made. If 
these relations are used, the following equations can be derived: m 0 = mp 0 (t 0 /t) 3 and o; 0 = 
Qj(pt/pot 0 ). By the use of these equations to eliminate m 0 and a 0 from Eq. 4.3, the following 
equation is obtained: 

log [a(tp/t 0 Po)] = J + B log [m(t 0 /t) 3 (p 0 /p)] ( 4 - 4 ) 

Since m = tps and a = sCp/m, Eq. 4.4 can be reduced to 

log (C D /t„p 0 ) = J + B log [potjftsVt) 2 ] (4-5) 

Since t 0 and p 0 were assumed to have constant values in the analysis expressed by Eq. 4.3, 

Eq. 4.5 can be simplified to 

log C D = K, + K 2 log (s fyt) (4.6) 

where K ( and K 2 are constants. Thus log C D is a linear function of log s Vt if log a can be 
assumed to be a linear function of log m (Eq. 4.3). 
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To verify the relation derived above, a plot was made of vs. s 1 / 2 /t, using all the data 
for plates obtained in this study. Figure 4.21 shows this plot, along with the results of a least- 
squares regression analysis. The data presented here are for 40 regular plates, referred to 
in Fig. 4.13, and 77 irregular plates, referred to in line segments 4 and 6 in Fig. 4.20. It is 
instructive to note that the drag coefficients for regular plates tend to be slightly higher than 
those for the irregular ones. Nevertheless, the standard error of estimate in C D for all plates 
was only approximately 5.6 per cent, being between that for window-glass fragments (3.4 per 
cent) and that for plate-glass fragments (7.8 per cent). The areas for both types of glass 
fragments were computed from mass, thickness, and the measured glass density (39.657 g/in. 3 ). 

4.5 SMALL-ANIMAL DATA 

Six mice, six rats, six guinea pigs, and five rabbits were dropped during these tests. The 
animals, which had been given lethal doses of Nembutal before the drops, were dropped as 
many as five times each. 

The animals were dropped in two different initial orientations: first, with one side toward 
the direction of fall (lateral) and, second, with the abdomen toward the direction of fall (ven¬ 
tral). Figure 4.22 is a photograph of the presented areas of a few of the animals in these 
positions. The areas of most of the animals were measured from similar photographs and are 
recorded in Fig. 4.23, with circles representing the lateral areas and squares representing 
the ventral areas. The straight line on this plot has a slope of 2 / 3 . It tends to show that, at 
least for the mice, rats, and guinea pigs, the area is proportional to the mass to the 2 / 3 power 
(see Ref. 5). It should be remembered that, although the animals tended to hold their initial 
orientations as they fell, the areas on Fig. 4.23 are larger than those “seen” by the wind since 
the hair tends to lie close to the body as the animal falls. 

Figure 4.24 is a plot of the experimentally determined a as a function of mass; the circles 
represent the a of animals dropped laterally and the squares represent the a of animals 
dropped ventrally. The regression line was computed using all the points. The slope of-0.3240 
is near the theoretical slope of — which was determined by assuming a constant density and 
a constant drag coefficient for any group of similar objects (see Secs. 4.3.3 and 4.3.4). 

Figure 2.3 shows the impact velocity as a function of drop height and a. The a values, 
determined experimentally from the animal data, are plotted on the theoretical line for the 
actual height of drop. In addition, the means and standard deviations of the experimentally de¬ 
termined a are tabulated in Fig. 4.24. They are shown in Fig. 2.3 and can be used to estimate 
the means and standard deviations of the impact velocities of the various animals. 


4.6 DRAG PROPERTIES FOR MAN 

The least-squares regression line in Fig. 4.24 can be extrapolated up to a mass of 168 lb 
(75,977 g) to give an estimation of the a for man (0.0269 ft 2 /lb). 

The extrapolated value for the a for man was compared to Schmitt’s report on wind- 
tunnel investigations of air loads on human beings 6 whose average weight was 168 lb. The data 
from this report were used to compute and average the a for a group of eight nude men in the 
lateral and ventral positions. A value of 0.0336 ft 2 /lb was obtained. This value is higher than 
the 0.0269 ft 2 /lb obtained from the extrapolation of the animal data, but this can be explained if 
the fact is taken into consideration that the animals tested were much more spherically shaped 
than man and thus have a lower a for the maximum presented-area position. 

Another way to estimate the a of man is by the concept of the total surface area. For a 
sphere the total surface area is 47 rr 2 and the average presented area is 7 rr 2 . For a cube the 
total surface area is six times the area of one face, and the average presented area of a cube 
in random orientation is 3 /> times the area of one face . 7 The ratio of the average presented 
area to the total surface area is ( 7 rr 2 / 47 rr 2 ) sphere = V 4 and [( 3 / 2 )/ 6 ] cube = 

In general, if an object has no cavities (ellipsoid, cube, etc.), its average presented area 
in random orientation is one-fourth the total surface area (see Appendix, Sec. A.3, for proof). 
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Since a man does have concavities, it is not actually the total surface area that should be con¬ 
sidered, but rather an effective surface area. 

Bohnenkamp’s work with effective radiating surface area of an erect man has shown this 
area to be from 83 to 85 per cent of the total surface area. 5 The total surface area can be ob¬ 
tained from various sources, 8-9 but for this study the figures given in the Schmitt report 6 were 
used (an average of 20.7525 ft 2 total surface area and 167.5 lb weight). With these data the 
average value of s/m was found to be x 84% x (20.7525 ft 2 /167.5 lb) = 0.0260 ft 2 /lb. If Cj) is 
assumed to be equal to unity, this figure (0.0260 ft 2 /lb) can be compared to the extrapolated 
value for a of 0.0269 ft 2 /lb (Fig. 4.24). 

The a values derived from the Schmitt report for various orientations of the men relative 
to the wind were used to calculate the average a for a rigid, clothed man, giving a value of 
0.0281 ft 2 /lb (see Appendix, Sec. A.4). The average drag coefficient for a rigid, clothed man 
would then be the average a divided by the average s/m, or 0.0281/0.0260 = 1.08. This value 
compares very well with the approximate values listed by Hoerner, 10 which include drag coeffi¬ 
cients of 1.0 for a parachutist in free fall, 1.0 for a ski runner on a slope, and 1.2 for a ski 
jumper in air. It should be remembered that the average a of man would be lowered if those 
positions where a man is doubled up, instead of being rigid, were also considered. 
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TABLE 4.1 — ALPHA VARIATION WITH MASS DERIVED FROM EXPERIMENTAL 

DATA FOR SPHERES AND STONES* 
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vs. sphere mass. The least-squares regression lines were computed from test data assuming constant Cq. 
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Fig. 4.5—Experimentally determined a vs. mass for painted gravel. Included are two least-squares regression lines prepared for two sets of 
conditions as stated on figure. 





46 


Fig. 4.6—Experimentally determined cc vs. mass for limestone fragments. Results of 
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from 10P4 and lSla locations at NTS are shown. 
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Fig. 4.8 Experimentally determined ct vs. mass for IS la natural stones. Included are results of the least-squares regression analysis for 
natural stones from 10P4 and IS la locations at the NTS. 
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Fig. 4.9—Experimentally determined cc vs. mass for natural stones from NTS at 4.3GTS location. 
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Fig. 4.11—A summary of the stone analysis showing the least-squares regression lines describing oc vs. mass for the various stone groups 
tested. 
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mass for military debris showing results of least-squares analysis. The theoretical line for 







Fig. 4.14 Drag coefficients of cubes (dropped with one face down) computed from their times of fall. 








Fig. 4.16 Front views of some of the larger pieces of window glass (0.125 in. thick) used in drop-tests. 
Numbers indicate masses in grams. 
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Fig. 4.17—Grouped data for window glass dropped flat showing means and standard deviations and least-squares regression lines describing 
oc vs. mass. This figure indicates the need for four straight line segments in the analysis. 





thick Window Glass Fragments 
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All Glass 



Section 1 10<m<220, log a = .55659 - - log m Section 3 I270<m < 1574, log a =-3.02323 + .89295 log m 

n = t54 cr = .09014 log units—20.9% n =10 cr = .04525 log units=10.4 % 

Section 2 220<m<1270, log a = .13874 - .036374 log m Section 4 1574<m<10, log a = -.28024 + .034985 log m 

n = 81 a - .06348 log units = 14.6 % n = 46 cr = .01473 log units= : 3.4% 

Edge Window 

3 




0.3 



Flat Plate 



10 10 2 10 3 10 4 10 s 
Section 1 10<m<860, log a = 55659 - -y log m Section 6 860<m< 3(10), log = - .49172 + .023921 log m 

n =155 O' = .09014 log units=20.9 % n = 31 cr = .03416 log units=7.8 % 


Fig. 4.20 —Summary of the glass-fragment data. This shows regression lines describing oc vs. mass for 
(1) all the glass, (2) window glass dropped flat, (3) window glass dropped on edge, and (4) plate glass 
dropped flat. 
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LATERAL POSITION 


VENTRAL POSITION 



MOUSE 22.45 g 





RABBIT 1829 g 


Fig. 4 ,22 —“Silhouette photographs of some of the animals used in the drop tests. The areas of the animals in the 
lateral and ventral positions were determined from full-scale photographs of this type. 
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Fig. 4.23—Silhouette areas vs. mass of animals used in the drop tests. These areas were photographically determined in the lateral and 
ventral positions. The line, which has a slope of 2 / 3 , is given for reference. 





log a = ,01153 — 32400 log m (gms) 







Chapter 5 

SUMMARY 


The purpose of this study was to determine certain aerodynamic-drag parameters of small 
irregular objects, such as window-glass fragments, stones, steel fragments, and small labora¬ 
tory animals. The drag parameters of interest were those necessary for the computation of the 
velocity vs. time history of the object when exposed to the winds associated with nuclear- 
produced blast waves. 1 These parameters were combined into one quantity called the accelera¬ 
tion coefficient (a) and defined as the product of the area presented to the wind and the drag 
coefficient divided by the mass. 

The experimental method used to determine a was to measure the time required for a test 
object to free fall a given distance. To relate a to drop time, it was necessary to assume that, 
first, the drag coefficient and, second, the area normal to the object's motion were both con¬ 
stant throughout the fall; i.e., a was constant. Experimentally, the second condition was not 
difficult to satisfy since this only required that the object not rotate. It is well known that the 
drag coefficient varies with wind speed; however, it was shown from the experimental results 
that the error introduced in the measurement of a by assuming that the drag coefficient was 
constant was not significant, at least in so far as irregular objects were concerned. The follow¬ 
ing derived relation was used to determine a from measured drop time: 

T = (2H/g) ,/2 (paH.)- 1 ' 2 cosh -1 e paH/2 
where T = drop time 

g = acceleration of gravity (assumed constant) 

H = height of drop 

p = air density (assumed constant) 

a = acceleration coefficient (assumed constant) 

The equipment used in the experiment consisted of a mechanism that would simultaneously 
drop the test object and start an electronic timer, an impact plate, and a contact microphone 
attached to the lower side of the impact plate. When the test object made contact with the plate, 
an electrical signal from the contact microphone stopped the timer. 

The dropping technique as used in these experiments, was explored by dropping several 
types of regular objects (spheres, cubes, and plates); the results obtained were compared with 
those obtained using other techniques. The results are presented in this report for irregular 
objects. The drop times for two types of spheres were also used to check the experimental 
equipment for consistency of operation at various times during the drop tests. 

Several types of stones were used as test objects. Most of the stones had been used at the 
NTS to study the production of secondary missiles by blast winds produced by nuclear explo¬ 
sions. 2,3 From 29 to 45 stones of each type were dropped, and a was determined for each. The 
results were analyzed by regression analyses to determine a for each type of stone, and the 
results are listed below. The units for a are ft 2 /lb and for mass (m), mg. Sigma (a) is used to 
designate the standard error of estimate in log units. 
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Gravel (painted): 


log a = 0.57610 - 0.38803 log m 
30 < m < 1500 a = 0.07463 

Limestone fragments: 

log a = 0.49926 - 0.32159 log m 
130 < m < 15,000 a = 0.08811 

Natural stones 3 from NTS at location 10P4: 

log a = 0.48759 - 0.32835 log m 
11 < m < 4100 a = 0.07063 

Natural stones 3 from NTS at location ISla: 

log a = 0.45171 - 0.31263 log m 

17 < m < 7100 a = 0.08684 

Natural stones 3 from NTS at location 4.3 GTS: 

log a = 0.73833 - 0.30758 log m 
11 < m < 2300 a = 0.1294 

Pumice fragments: 

log a = 0.63748 - 0.25466 log m 

18 < m < 2300 cr = 0.08948 

Fragments of steel (mostly shrapnel) were also used as test objects, and the results were 
analyzed in a manner similar to that described for stones. The resulting regression equation 
was 

log a =-0.74223 - 0.11315 log (m/1000) 

4500 < m < 150,000 a = 0.09974 

Fragments of ordinary window glass 0.125 in. thick were dropped, both with a plane sur¬ 
face down and edge-on, to determine the maximum and minimum a as a function of mass. In 
addition, plate-glass fragments 0.225 in. thick were dropped with the plane surface down. For 
small fragments the a values determined were independent of both orientation and type of glass 
(ordinary window glass or plate glass). It was found that the relation between a and mass 
could be represented as a straight line on logarithmic paper, but only for definite mass in¬ 
tervals. The regression equations derived from the experimental data are given below (n 
represents the number of fragments dropped; the remaining nomenclature is the same as that 
used previously). 

Window glass and plate glass dropped flat or edgewise for 10 < m < 220: 

log a = 0.55659 - (1/3) log m 
n = 155 a = 0.09014 
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Window glass dropped flat for 220 < m < 1270: 

log a =-0.13874 - 0.036374 log m 
n = 81 a = 0.06348 

Window glass dropped flat for 1270 < m < 1574: 

log a = —3.02323 + 0.89295 log m 
n = 10 a = 0.04525 

Window glass dropped flat for 1574 < m < 10 5 : 

log a = -0.28024 + 0.34985 log m 
n = 46 a =0.01473 

Window glass dropped edgewise for 220 < m < 5000: 

log a = 0.12689 - 0.14982 log m 
n = 88 a = 0.1013 

Plate glass dropped flat for 860 < m < 3 x 10 5 : 

log a = 0.49172 + 0.023921 log m 
n = 31 a = 0.03416 

A useful formula was derived which related drag coefficient (Cd) for plates of random 
shapes and finite thickness to the ratio of the square root of plate surface area (s) to thickness 
(t). The drop-test data used in this analysis included that for window- and plate-glass frag¬ 
ments dropped flat as well as that for plates of regular shapes (square, circular, triangular) 
made of brass and plastic. The regression equation noted below was computed from data for 
117 test objects whose shape factor (s 1/2 /t) varied from 2 to 65. 

log C D = 0.019445 + 0.048498 log (s V2 /t) 
a = 0.02428 

(Note: s 1/2 and t are expressed in the same units.) 

Shadowgraphs were prepared for mice, rats, guinea pigs, and rabbits in ventral and lateral 
positions to determine their presenting areas. In general, the areas determined (Fig. 4.23) 
were proportional to the mass of animal raised to the 2/3 power. The animals were then 
dropped ventrally and laterally to determine their a. It was found that there was no significant 
difference in a determined for these positions. The experimentally determined a values for 
each animal species were: 

Mice (20.44 g) 0.3769 ± 0.0425 ft 2 /lb 

Rats (184.48 g) 0.1853 ± 0.0142 ft 2 /lb 

Guinea pigs (530.1 g) 0.1547 ± 0.0196 ft 2 /lb 
Rabbits (2087 g) 0.0794 ± 0.0155 ft 2 /lb 

These data were used to determine a as a function of animal mass. 

log a = 0.01153 - 0.32400 log (m/1000) 
n = 84 <j = 0.06496 


An estimate of the a for a 168-lb man was determined to be 0.0269 ft 2 /lb by the above 
equation. Another estimate was made for the a of a 168-lb man using the results of a wind- 
tunnel study reported by Schmitt. 4 For positions corresponding to those of the dropped animals 
(ventral and lateral), an average a of 0.0336 ft 2 /lb was obtained. 

A method for estimating the average presenting area of any object in random orientations 
was developed. It was shown that the average presenting area is one-fourth the total surface 
area, provided no cavities exist. To approximate the no-cavity condition, the radiating surface 
area for man was used to estimate the average a for the human subjects used in the wind- 
tunnel study mentioned above. Assuming a drag coefficient of 1.0, a value of 0.0260 ft /lb for a 
was obtained. This value compared well with a value of 0.0281 ft 2 /lb obtained by averaging the 
wind-tunnel data for men in various positions relative to the wind. 


REFERENCES 

1. I. Gerald Bowen et al., A Model Designed to Predict the Motion of Objects Translated by 
Classical Blast Waves, Report CEX-58.9, June 1961. 

2. I. Gerald Bowen, Allen F. Strehler, and Mead B. Wetherbe, Distribution and Density of 
Missiles from Nuclear Explosions, Project 33.4, Operation Teapot Report, WT-1168, 
December 1956. 

3. I. Gerald Bowen et al., Secondary Missiles Generated by Nuclear-produced Blast Waves, 
Project 33.2, Operation Plumbbob Report, WT-1468, in preparation. 

4. Thomas J. Schmitt, Wind-tunnel Investigation of Air Loads on Human Beings, Report DTMB- 
892, The David Taylor Model Basin Aerodynamic Laboratory, January 1954. 


68 


Appendix 

DERIVATION OF EQUATIONS 


A.l EQUATIONS OF MOTION 

A. 1.1 Integration of the Basic Equation 
Equation 2.9 can be integrated 


to give 



(A.l) 


(A.2) 


Now, solving Eq. A.2 for v 2 gives 


v 2 = 



(» 



(A. 3) 


where h = h 2 - hq represents the height fallen from the point where the velocity equalled v t . 
Substituting v 2 = dh/dt in Eq. A. 3 and integrating from zero height of drop to h 



gives 


*2 “ ti 


(f)“ 
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Substituting h = (h 2 - h t ) from Eq. A.2 into Eq. A.5 gives 



(A.4) 


(A.5) 


(A. 6) 
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A. 1.2 Dimensionless Analysis 


Equations 2.4 and 2.8 of the text can be written in the form 



(A.7) 


where f(vD/z^) = f(Rd) = C^j, the drag coefficient. 

Now, g, D, and v were considered constant so that they could be placed inside the differ¬ 
entials. Algebraic manipulations were then made on Eq. A.7 to make dimensionless quantities 
as follows: 



Defining Y 
Z 
K 


= hg(D/V ) 2 (distance numeric) 

= vD/v (velocity numeric) 

= (p/g) (s/m) (iv/D) 2 (acceleration-coefficient numeric) 


(A. 7a) 


as numeric parameters, the above dimensionless equation becomes 



= 1 - |kz 2 f(Z) 


(A.8) 


By basic definition, V = AH/AT or AT = AH/V, Eq. A.8 does not contain the time numeric. 
However, the incremental step of time numeric, defined as AX, would be the incremental step 
of distance numeric divided by the velocity number. Restated in terms of an equation, this is 


Thus 


AX - 


AY Ah g ( D/u ) 2 D Ah D 

T s v(DA) =g V~" g V AT 


Y-n- D t 

X - g t 
v 


(time numeric) 


(A.9) 


Any previously derived equation can now be converted to the nondimensional form by re¬ 
placing the primary parameters v, h, t, and s/m, by the numerics Z, Y, X, and K. The p, v, 
and g quantities are dropped in the nondimensional form as they go into the numerics. Cp re¬ 
mains the same in the nondimensional equations. 

Equations A.2 and A.6 become 
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(A.10) 


and 
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(2KC d ) 1/s 


In 



z t l 

(KCrf 

l 2 i 

72-t 

1] 


1 + Z 2 ' 



Z 2 I 

i'KCd') 
< 2 } 

H 


1 + z,l 

^kcdV^ 
2 / . 

" > 


(A.11) 
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(T — T ND ) or time-of-fall minus no-drag time-of-fall parameter* (Sec. 2.5), would now be 

(A.12) 


T - T nd = T - ~ X — (2y/ 2 


Another useful parameter (Sec. 2.5) would be 


T - T 


ND 


--i2"T ND -T \ ) V + 12 


X 


-(2Y) I/2 (l +-§KY) 


(A.13) 


If the approximate average value of the drag coefficient of an object is substituted for C d 
in Eq. A.13, the parameter then represents the scaled time difference between the actual time 
above no-drag time and an approximation to the time above no-drag time considering the drag 
coefficient to be constant at C]> It should be noted one can go from the variable of Eq. A. 13 to 
that of Eq. A. 12 by adding 2Y(C D /12) KY, which is the approximation to the time of fall minus 

no-drag time of fall in nondime ns ional form. 

The no-drag equations, Eqs. 2.18 and 2.19, can likewise be put in numeric form to give 


Z 

Y 


ND “ X ND 

_ (Xnd) 2 

ND-2 


(A.14) 
(A. 15) 


Equations A. 10 and A. 11 are the equations actually used to predict the accurate drop times 
of two spheres (Sec. 4.2). The C D curve was broken up into about one hundred steps over the 
Reynolds-number range encountered by these spheres, and the Cj} values were used to integrate 
step by step to the total height of drop. 


A.2 TERMINAL VELOCITY 


Although not apropos of the present experiment, it should be noted that a useful relation 
can be derived from Eq. 2,13. As a body falls farther and farther, its velocity tends to reach a 
finite value called terminal velocity. This terminal velocity can be found from Eq. 2.13 by 


lim V = lim (2gH)^ 2 

H-*°° H -* 00 



v 2 



(A.16) 


Figure A.l is a plot of terminal velocity vs. a for two different conditions of p and g. 

The percentage of terminal velocity reached by an object with constant ct after it has fallen 
a distance H can be computed from Eqs. 3.13 and A. 16 by 


%V 


terminal 


100 


V actual 
^terminal 


100 


y 2 

(2gH) /2 


/1 _ e -paH\ 1/2 
\ pa H / 



= ioo (i - e - p “ H ) 1/2 


(A.17) 


Figure A.2 is a plot of the above relation showing the percentage of terminal velocity as a 
function of a and H. Equation A. 16 is independent of any initial conditions and represents the 


*The subscript ND is used to indicate the no-drag value of a parameter; i.e., T ND is the 
time an object would take to fall, assuming the object did not experience drag force. 
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velocity to which the instantaneous velocity is tending at any given time; whereas, Eq. A.17 
requires that the velocity be zero when the height is zero and also, that p, a, and g remain 
constant throughout the fall. 

A.3 AVERAGE PRESENTED AREA 

Consider a smooth surface, S, which has no cavities (part a, Fig. A.3); i.e., a plane that is 
tangent to S at any point does not cut through the surface. Then at some point on S construct a 
tangent plane that approximates S over a small area dA in the vicinity of the point of tangency. 
Let N be the normal to the plane at the point of tangency, E be oriented in any reference direc¬ 
tion fixed in space, and 9 be the angle between N and E. As N takes random orientations, not 
all values of 9 are equally likely. To find the probability distribution of 0, one compares the 
area on the surface of a sphere between 9 and 9 + d0 to the total area of the sphere. That is to 
say, consider a sphere whose center is always located at the point of tangency of the plane on 
the smooth surfaces. Now N can intersect the sphere in any point on the sphere, and all such 
points are equally likely. Then the probability that N will intersect the sphere in a specified 
area of the sphere is just the ratio of this area to the total area of the sphere. In the case 
under consideration, the specified area is the total area on the sphere between the angular 
limits of 9 and 9 +'d0 with respect to the fixed reference direction E. This area (see part b, 
Fig. A.3) is 27r(r sin 9) r d 9; so the probability that the angle between N and E is between 9 and 
9 + d0 is 


P(0) d0 - 


2?r(r sin 9) r d0 
47rr 2 


sin 9 
2 


d 9 


(A.18) 


The presented area in the direction of E produced by dA (part a, Fig. A.3) is 


Presented area = dA cos 9 


when 0 < 9 < 


7 T 


= 0 


when ~< 9 < tt 


(A.19) 


For 7t/2 < 9 < tt , dA could not be seen from the direction of E because it is behind S. The aver¬ 
age presented area would then be the integral of the presented area as a function of 9 weighted 
with the probability distribution of 9, or 


Average presented area = 



sin 9 
~2~ 


dA cos 9 d 6 


sin 2 9 
4 


dA 


tt/ 2 = dA 
o 4 


(A. 20) 


For the whole surface S, the average presented area divided by the total surface area would be 


Average presented area 
Total surface area 


fsT L _lfs dA _l 


(A.21) 


A.4 AVERAGE a FOR A RIGID MAN 


Two angles were used to fix the man with respect to the wind, 9 and cp (part c, Fig. A.3). 
Alpha can be computed directly from Schmitt’s report 1 for 9 = 90°; 0° < (p < 360° (upright 
position) or for 0° < 9 < 180°; <p = 90° (side to the wind). For intermediate values the following 
formula was used to approximate ol 


a(9,(p) = a(0°) + 


[(*(90°,0) - «(0 0 )] [a(9, 90°) - of(0°)] 
01(90°,90°) - a{0°) 


(A.22) 
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where Q{(0 o ) represents the cn of a man head-on to the wind (0 — 0) and is the same irrespective 
of <p , similarly q!( 180 o ) represents the a of a man feet-on to the wind (0 = 180°) and also is 
constant for all values of fa From Schmitt’s report 1 0(0°) = 0(180°) for a clothed standing man 
and has a value of 0.00627 ft 2 /lb. Thus Eq. A.22 holds exactly for all the conditions that can be 
computed directly from Schmitt’s report; i.e., substituting in Eq. A.22 the various values of 0 
and <p for which a is known, the following relations for a are obtained: 


a(0°,<p) = q?( 0°) 
a(180°,<p) = a(0°) 
o:(90 ° f (p) = a(90°,<p) 
a(0,9O°) = a(0,9O°) 


(A.23) 


Equation A.22 was arrived at by assuming [a(Q,fa) - a(Q)]/[a(6, fa) - a(0)] is constant for all 
values of 0 when fa and fa are fixed or, alternatively, by assuming [ot(6 u <p) - a( 0°)]/[a{e 2 ,fa - 
a(0°)] is constant for all values of <p when 9\ and 02 are fixed. In connection with the first 

assumption, one can write 


oi(6,<p) - a{0°) _ a(9 0°,<ft) - o?(0°) (A 24) 

o:(0,9O o ) - a(0°) a(90 o ,90 o ) - qi( 0 o ) 

for all values of 0 and <p. Solving this for a(0,<£) gives Eq. A.22. The second assumption will 
yield Eq. A.22 in a similar manner. Since 0 and (p are independent variables, the average a 

would be given by 


a(6,<p) = o?(0°) + 


[a(90°,(p) - q?(O°)1[q!( 0,9 O o ) - a(0°)] 
a(90 o ,90°) - a((f) ~ 


(A.25) 


where the bars indicate average values. 

The value q!( 0,9O 0 ) was computed by weighing ckf(0,9O°) with the probability distribution 
(sin 0/2) d0 (Eq. A.3) and integrating graphically from 0 = 0 to 0 = n radians. The value for 
a(90°,<p) is just the average value of a(90°,(pj with no weighing (all values of <p are equally 
likely), and a(90°,90°) and a(0°) are particular a values computed directly from Schmitt’s re¬ 
port. The value for $(0,(£>) computed from Eq. A.25 represents the average oi for a rigid, 
straight man, assuming every possible orientation with respect to the wind is equally probable. 
(Average a computed from Schmitt’s report is 0.0281 ft 2 /lb.) 


REFERENCE 

1. Thomas J. Schmitt, Wind-tunnel Investigation of Air Loads on Human Beings, Report DTMB- 
892, The David Taylor Model Basin Aerodynamic Laboratory, January 1954. 
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700 



Terminal velocity vs. a for standard conditions at sea level (curve I) and typical conditions at 5300-ft altitude (curve II) 
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Fig. A.2— Percentage of terminal velocity vs. height of fall and a. When p = 0.0653 lb/ft 3 (typical air density for 5300 ft above sea level), 
the height can be read directly, if p does not have this value; the height variable must be multiplied by the appropriate factor to convert to 

actual height of fall. 
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CIVIL EFFECTS TEST OPERATIONS REPORT SERIES (CEX) 

Through its Division of Biology and Medicine and Civil Effects Test Opera¬ 
tions Office, the Atomic Energy Commission conducts certain technical tests, 
exercises, surveys, and research directed primarily toward practical applica¬ 
tions of nuclear effects information and toward encouraging better technical, 
professional, and public understanding and utilization of the vast body of facts 
useful in the design of countermeasures against weapons effects. The activities 
carried out in these studies do not require nuclear detonations. 

A complete listing of all the studies now underway is impossible in the 
space available here. However, the following is a list of all reports available 
from studies that have been completed. All reports listed are available from 
the Office of Technical Services, Department of Commerce, Washington 25, 

D. C., at the prices indicated. 

CEX-57.1 The Radiological Assessment and Recovery of Contaminated 
($0.75) Areas, Carl F. Miller, September 1960. 

CEX-58.1 Experimental Evaluation of the Radiation Protection Afforded by 
($2.75) Residential Structures Against Distributed Sources, J. A. Auxier, 

J. O. Buchanan, C. Eisenhauer, and H. E. Menker, January 1959. 

CEX-58.2 The Scattering of Thermal Radiation into Open Underground 

($0.75) Shelters, T. P. Davis, N. D. Miller, T. S. Ely, J. A. Basso, and 

H. E. Pearse, October 1959. 

CEX-58.7 AEC Group Shelter, AEC Facilities Division, Holmes & Narver, 
($0.50) Inc*, June 1960. 

CEX-58.8 Comparative Nuclear Effects of Biomedical Interest, Clayton S. 
($1.00) White, I. Gerald Bowen, Donald R. Richmond, and Robert L. 
Corsbie, January 1961. 

CEX-58.9 A Model Designed to Predict the Motion of Objects Translated by 
($1.25) Classical Blast Waves, I. Gerald Bowen, Ray W. Albright, E. Royce 
Fletcher, and Clayton S. White, June 1961. 

CEX-59.1 An Experimental Evaluation of the Radiation Protection Afforded 
($0.60) by a Large Modern Concrete Office Building, J. F. Batter, Jr., 

A. L. Kaplan, and E. T. Clarke, January 1960. 

CEX-59.4 Aerial Radiological Monitoring System. I. Theoretical Analysis, 
($1.25) Design, and Operation of a Revised System, R. F. Merian, 

‘ j t g. Lackey, and J. E. Hand, February 1961. 

CEX-59.13 Experimental Evaluation of the Radiation Protection Afforded by 
($0.50) Typical Oak Ridge Homes Against Distributed Sources, T. D. 
Strickler and J. A. Auxier, April 1960. 


